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Abstract. Two dense dark cloud cores representing dif- 
ferent stages of dynamical evolution were observed in a 
number of molecular spectral lines. One of the cores, Cha- 
MMSl in the Chamaeleon cloud i contains a Class proto- 
star, whereas the other, CrA C in the R Coronae Australis 
cloud, is pre-stellar. The molecules selected for this study 
are supposed to show significant abundance variations in 
the course of the chemical evolution. 

We find that the cores have very different chemi- 
cal compositions. Cha-MMSl exhibits characteristics of 
so-called 'early-type' chemistry with high abundances of 
carbon-chain molecules such as HC3N, CH3CCH and c- 
C3H2. However, it also has a large N2H"'" abundance, 
which is expected only to build up at later stages. In con- 
trast, none of the carbon-chain molecules were detected 
in CrA C. On the other hand, CrA C has a higher SO 
abundance than Cha-MMSl, which according to chem- 
istry models implies that it is chemically 'older' than Cha- 
MMSl. The most striking difference between the two cores 
is seen in the HC3N/SO abundance ratio, which is at least 
three orders of magnitude higher in Cha-MMSl than in 
CrA C. This result is somewhat surprising since starless 
cores are usually thought to be chemically younger than 
star-forming cores. 

Because of the high N2H+ abundance, we suggest that 
Cha-MMSl represents the 'late-time cyanopolyyne peak' 
that is predicted to occur when heavy molecules start to 



freeze onto grain surfaces (Ruffle et al. 1997). This would 
also be a more natural explanation for the carbon-chain 
molecules than the 'early-time' picture in view of the fact 
that the core is presently collapsing to form a star. The 
abundances observed in CrA C can be explained either 
by pure gas-phase models at late stages of evolution, or 
by the 'SO peak' which follows the second cyanopolyyne 



peak (Ruffle et al. 1999). Thus, the dynamical evolution 
in CrA C seems to have been very slow compared with 
that of Cha-MMSl, and we discuss possible reasons for 
this. 

We detected two SO emission maxima around Cha- 
MMSl, which lie symmetrically on both sides of the core. 



approximately on the line connecting the centre of Cha- 
MMSl and the position of Hcrbig-Haro object HH49/50. 
These SO peaks may signify the lobes of a bipolar out- 
fiow, and the obser vation supports the suggestion by 
Rcipurth et al. 1996 that Cha-MMSl is the central source 



of HH49/50. 

Key words: molecular processes - ISM: abundances, 
clouds, molecules - individual: CrA C, Cha-MMSl 



1. Introduction 

To understand the first steps of the star-formation pro- 
cess, we must be able to follow the evolution of molecu- 
lar cloud cores. One approach is to conduct spectral-line 
observations of molecular species whose abundances vary 
significantly with time. Conversely, combining such obser- 
vations with other types of diagnostics of cloud evolution, 
the chemical models may be tested. 

Examples of species that are believed to be useful for 
the aforementioned purposes are complex hydrocarbons 
and cyanopolyynes0. Their abundances are expected to 
peak at early stages o f the evolution of dark clouds (e.g. 
Herbst fc Leung 1989|) . This is because their production 
requires a high abundance of free neutral carbon (C i) ; 
a condition that is satisfied at early times in the chemi- 
cal evolution, whereas at later stages practically all free 
carbon has been exhausted in the production of CO. In 
contrast, molecules such as NH3 and SO increase in abun- 
dance with time, because their formation involves rela- 



tively slow neutral- n eutral reactions ( Suzuki et al. 1992 



Hirahara et al. 1995 ; Nilsson et al. 200C ), or because they 
(e.g. SO) are destroyed by Ci or C+. 

The validity of these model predictions have, to date, 
been tested in a rather limited sample of dark cloud 
cores. Most observations have been confined to the Tau- 
rus Molecular Cloud 1 ( TMC-1), L134A, and to a num- 



ber of "Myers' cores" ( [Pratap et al. 1997t |Swade 1989 
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Myers & Benson 1983, Benson et al. 1998; for a review 
see Ivan Dishoeck fc Blake 199S| ). However, for current 



chemistry models to be rigorously tested, the cloud sample 
needs to be more representative of the ISM. In the present 
paper we extend the sample of well-studied objects by two 
potentially very young molecular cloud cores. 

The first core selected for this study is located in 
the region of the reflectio n nebula Cederblad 110 in the 
Chamaeleon dark cloud i. Rcipurth et al. 1996 detected a 
strong 1.3mm dust continuum source towards this core, 
and designated it as Cha-MMSl. The core belongs to an 
extensive, northeast-southwest directed filament of dense 



gas near the centre of this cloud (Mattila et al. 198£). The 
second core, CrA C ( Harju et al. 1993| ) is located in the 
southeastern 'tail' of the R Coronae Australis molecular 
cloud. The distances to the Chamaeleon I and the R Coro- 
nae Australis clouds are 150 and 170 pc, respectively 
(|Knude fc H0g 19981). The smallest diameters of the C^^O 



cores are 0.11 and 0.05 pc for CrA C and Cha-MMSl, re- 
spectively. The C^^O column densities towards both cores 
are similar: 7V(C^®0) ~ 210^^ cm~^, which corresponds 
to a molecular hydrogen column density of A^(Il2) ~ 
10^2 cm-2 dHarju et al. 1993| ; pVIattila et al. 1989| ). The 
average H2 densities derived from these numbers assum- 
ing prolate geometries are ~ 3 10^ and ^ 7 10^ cm"*^ for 
CrA C and Cha-MMSl, respe ctively. The kineti c temper- 
ature in CrA C is 10 K or less ( Harju et al. 1993 ), whereas 
in Cha-MMSl the gas kinetic temperature as derived from 
ammonia observations is slightly elevated with respect to 
the typical dark cloud value, being ~ 14 K (Toth et al. 
2000, in preparation). 

Even though both cores are supposed to be young, 
they represent different stages of cloud evolution. 
Cha-MMSl probably contains a Class protostar 



(Reipurth et al. 1996), which has been recently iden- 
tified with a far-infrared source, designated CedllO 
IRS 10, by [Lehtinen et al. 2000| . The core lies next 
to a small group of newly-born stars, the closest 
being Ced 110 IRS 4 ( iPrusti ct al. 19911) and the 
newly detected IRS 11 ( Lehtinen et al. 200C ), which 
belong to the infrared source Class I. The core 
is associated with a compact, north-south oriented 



molecular outflow ( 


Mattila et al. 1989 


Prusti et al. 1991 


Lehtinen et al. 20001). 


Reipurth ct al. 1996 suggested that 



Cha-MMSl contains the driving source of the outflow and 
the Herbig-Haro objects HH 49/50 located some 10' south 
of the core. 

CrA C shows no sign of star formation activity. The 
core appears to be surrounded by an extensive halo: it 
is located near the centre of a neutral atomic hydrogen 



ring with a diameter of 40' (Llewellyn et al. 1981) which 
may represent the outer boundary of a massive Hi to 
H2 conversion region. Observations of continuum emission 
at 100 ^m and 200 ^m wavelenghts with the ISOPHOT 
instrument on the Infrared Space Observatory (ISO) re- 
veal a 200/zm excess coinciding with the centre of CrA C 



(Heikkila et al. 2000, in preparation), suggesting that a 
very cold clump is located inside this core. 

Based on these facts, one could expect CrA C to ex- 
hibit the characteristics of 'early time' chemistry, e.g. high 
abundances of unsaturated carbon compounds, whereas 
Cha-MMSl could be thought to have reached a more 
evolved stage. In order to test these ideas, we have es- 
timated for these dense cores the column densities of sev- 
eral molecules which show significant time dependence in 
chemistry models. 

The selected molecules are discussed in the next Sec- 
tion. In Sect. 3 we describe the observational procedure. In 
Sect. 4 the results from the observations are presented and 
these are further discussed in Sect. 5. Finally, in Sect. 6 
we summarize our conclusions. 



2. Observed molecules 

We discuss here the spectroscopic and chemical prop- 
erties of the molecules included in this study. Table |^ 
lists their chemical formulae, structure types (hfs means 
that the lines have hyperfine structure; o/p and A/E 
that the molecule has ortho and para or ^-symmetry 
and i?-symmetry forms, respectively), permanent electric 
dipole moments (/^electric) along the relevant molecular 
axises, and rotational constants {A, B and C). The 
data were extracted m ainly from the JPL catalogue^ 
( Poyntcr fc Pickett 1985| ); exceptions are the CH3OH 
data ( Sastry ct al. 1981 ; Anderson et al. 1990|) and the 
value for the dipole moment of CCS ( Murukami 199C ). 



2.1. Spectroscopic properties 

Carbon monosulphide, CS, is a linear molecule with 
a simple rotational spectrum; the states are denoted 
by the total angular momentum quantum number J. 
The simplicity of the spectrum is due to the facts that 
in the electronic ground state, the electrons in the CS 
molecule have paired spins and possess no orbital angular 
momentum. Hence, the ground state is denoted ^S. In 
addition, the C and S nuclei have no spin. 

The dicarbon sulphide radical, CCS, and the sulphur 
monoxide radical, SO, have more complicated spectra 
than CS. The electrons in the electronic ground states 
of CCS and SO have no net orbital angular momentum, 
but the total spin is non-zero due to two electrons 
having unpaired spins. Thus, the spin quantum number 
is 5 = 1 and the electronic ground state is denoted •^E 
(e.g. pordy fc Cook 1970| , Sect. 4.2). The electron spin 
is coupled to the weak magnetic field arising from the 
rotation of the molecule. Since 5=1, for each molecular 



^ An updated versio n is available on the WWW at 
http://spec.jpl.nasa.gov 
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rotational level with the rotational angular momentum 
quantum number N > 0, the total angular momentum 
quantum number, J, has three possible values: J = iV — 1, 
J = N or J = N+1. Consequently, the rotational levels 
> arc split into triplets, whereas for A^=0 only one 



a single observation of a J — > J — 1, multiplet (the 
lines lie closely in frequency), is considered to be a 
good approximation of the kinetic temperature (e.g. 
Bergin et al. 1994|). 



state (Jat ~ Iq) exists (e.g. Yamamoto et al. 1990) 



Hydrogen isocyanide, HNC, is a linear molecule. The 
nitrogen nucleus has a non-zero spin, and therefore 
possesses an electric quadrupole moment. This interacts 
with the electric field gradient of the molecule, which 
depends on the rotation. Consequently, the nuclear spin 
(/ = 1) is coupled to the molecular rotational angular 
momentum, J, to yield the total angular momentum, F 
( iGordy fc Cook 1970| , Sect. 9.4). The resulting hyperfine 
splitting of the rotational lines of HNC is smaller than for 
its isomer HCN. For example, for the lowest rotational 
transition (J = 1 ^ 0) the F = 1 ^ 1 and F = ^ 1 
components are separated from the main component 
f = 2 ^ 1 by -0 .27 and +0.41 kms~^ respectively 
( Frcrking et al. 1979| ). Hence, the hyperfine structure 
is usually not resolved in observations due to Doppler 
broadening in the source and the insufficient spectral 
resolution of the spectrometers available. 



Cyanoacetylene, HC3N, is also a linear molecule. Due 
to its high moment of inertia, it has a small rotational 
constant (B) and its rotational transitions J ^ J — 1 lie 
relatively closely in frequency. Many of these lines are 
located in the commonly observed frequency bands, and it 
is often used in multi-transition studies. Like HNC, it has 
hyperfine structure due to the nuclear spin of nitrogen. 
For higher J- values, however, the relative energy differ- 
ences between different i^-levels are very small, and so 
no hyperfine structure is seen in the spectra observed here. 

Methyl acetylene, CH3CCH, is a low dipole moment, 
symmetric top molecule. Its energy levels are described by 
two quantum numbers J and K. The former represents 
the total angular momentum and the latter is the pro- 
jection of J on the symmetry axis of the molecule. The 
lowest energy state in each i^-ladder is J—K. CH3CCH 
exists in two chemically different forms, the so called A 
and E symmetry species, depending on the relative ori- 
entations of the spins of the three hydrogen nuclei in the 
CH3 group. For the A-symmetry species, the K quantum 
numbers are multiples of 3 (i.e. K — Zn] 7i=0,l,2,...), 
while the £^-symmetry species has K = 3 n + 1 and 
3n -f 2; n=0,l,2,... . The ground-state energies for A 
and E CH3CCH differ by 8.02 K, E lying higher in 
energy. There are no allowed electric dipole transitions 
between two states belonging to different iC-ladders. 
Consequently, the relative populations of two /^-ladders 
are determined by collisional transitions and therefore 
depend on the kinetic temperature. In fact, the so-called 
'rotational temperature', which can be determined from 



Cyclo-propenylidene, C-C3H2, is one of the few cyclic 
molecules detected in interstellar space. It is a slightly 
asymmetric, almost oblate top. C3H2 has ortho- and 
para-states originating from the two possible relative 
orientations of the spins of the two hydrogen nuclei 
(having ^ = 5 each). In the ortho-state the spins are 
parallel (/ — 1), whereas in the para-state the spins are 
anti-parallel (/ — 0). Since the ortho- and para- forms 
cannot be interchanged in electric dipole transitions or 
easily in chemical reactions, they can be seen as two 
different chemical compounds. The ground-state of the 
ortho-form {JKa,K^ = lo,i) lies 2.35 K higher than that 
of the para- form [Jk^.k^ = Oo,o)- The electric dipole 
transitions of C-C3H2 are of 6-type. This means that 
the quantum numbers describing the projection of J on 
the molecular a— axes and c— axes, Ka and Kc^ change 
according to the rule Ai^a — ±1, Ai^c — ±1 (see 
Vrtilek et al. 1987| ). 



Diazenylium, N2H+, is a linear molecular ion. Its 
rotational transitions have hyperfine structure due to 
the coupling of the spins of the two nitrogen nuclei 
{I1J2) to the rotational angular momentum (J). The 
energy levels are labelled using the quantum numbers 
J, Fi and F, where F i ^ J + h and F = Fi ^ h 
(e.g. ICaselli et al. 1995| ; pordy fc Cook 1970t Sect. 9.5). 
The J = 1 — > transition observed here has seven well 
resolved hyperfine components. 

Methanol, CH3OH, is a slightly asymmetric top (e.g. 
Townes & Schawlow 1975, chapters 4 & 12). It has a 
rather complicated rotational spectrum, since it exhibits 
so called hindered internal rotation: the hydroxyl (OH) 
group can rotate with respect to the methyl (CH3) group. 
This gives rise to two different forms of CH3OH between 
which both radiative and collisional transitions are for- 
bidden: the A-symmetry and F-symmetry forms. Chemi- 
cal reactions that might convert A to i? or vice versa are 
believed to occur only on very long time scales. There- 
fore the two forms of CH3OH can be treated as chemi- 
cally distinct species. The ground-state of the F-symmetry 
form (Jfe = l_i E) lies 6.95 K higher than that of the A- 
symmctry form [Jk = Oq A+). The relevant quantum num- 
bers describing the rotational spectrum are J and fc, with 
k = Ka — Kc- The electric dipole transitions observed here 
are of a-type, i.e. they obey the selection rule A/c ~ 0. Like 
other symmetric or slightly asymmetric tops, CH3OH may 
provide useful information about the gas temperature. 
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2.2. Chemistry 

The buildup of cyanopolyynes proceeds mainly via re- 
actions between complex hydrocarbon ions and nitrogen 
atoms (Herbst & Leung 1989). HC3N is the simplest 



cyanopolyyne, and its abundance correlates well with 
heavier cyanopolyynes (e.g. Federman et al. 199C). Due 
to its large dipole moment, as well as being chemically an 
early-time species, HC3N should be a useful probe of the 
conditions in pre-stellar cloud cores. 

The main reactions leading to HC3N and 
CH3CCH involve either the CciH:^ or C3H;]' ions 



(Huntress & Mitchell 1979). The neutral reactions 
involving C2H2 and CCH suggested recently by 



[Turner et al. 1998 and Turner et al. 1999 do not change 
this fact, since these species are also derivatives of the 
aforementioned molecular ions. Therefore, HC3N and 
CH3CCH abundances are expected to be well correlated. 



Reactions of nitrogen-bearing ions, such as HCNH+, 
with neutral hydrocarbons may also contribute to the 
production of cyanopolyynes ( Herbst fc Leung 19891 ). 
Since HCNH + is the precursor ion of HCN and HNC 
(Herbst 1978), a correlation between these isomeric 



molecules and cyanopolyynes is possible. In this paper, 
we shall not address the as yet unresolved problem 
of the variation of the HNC/HCN abundance ratio 
from source to source (for a recent overview see e.g. 
Talbi & Herbst 19981). We instead used HNC merely to 



survey the extension of the core CrA C, where HC3N was 
not detected. This should be as abundant as HCN but 
easier to detect due to the smaller hyperfine splitting. 



that the production of NH3 involves molecular nitrogen, 
N2, which is formed in neutral- neutral reactions and 
becomes abundant in the later stages of cloud evolu- 
tion when carbon chain molecules have already been 
lost in reactions with ions such as He"*", H+ and Hj 



(Suzuki et al. 1992). N2H"'" is also formed from N2 by 



N2 



H 



NoH^ 



Ho and its abundance has been 



found to follow closely that of NH3 (Hirahara et al. 1995 



Nejad fc Wagenblast 199£). Therefore, if the CCS/NH3 



abundance r atio can be used as an age indicator as 
suggested by Suzuki et al. 1992 , the same should be true 
for CCS/N2H+. The depiction of CO onto grain surfaces 
further increases the N2H"'" abundance by making H^ 
available for non-carbon bearing species to react with 
( [Nejad fc Wagenblast 1999| ). 



SO is destroyed by neutral carbon and, like N2H+, it 
is formed in neutral-neutral reactions and benefits from 
the freezing-out of CO. Its abundance therefore increases 
slowly and probably remains low as long as the gas is 
rich in neutral atomic carbon. The SO abundance first 
increases when carbon is consumed in the production 
of CO, i.e. at late s t ages in the chemical evo lution 
( [Bcrgin fc Langer 1997j ; [Nejad fc Wagenblast 1999[) . The 



late peaking of SO and the stability of CS (until it 
too is frozen out) has led to the suggestion that the 
SO/CS abundance could be used as a chemical clock 
( [Ruffle et al. 1999D . [Nilsson et al. 2000[ came to the same 
conclusion on the basis of pure gas phase chemistry 
models. However, the SO/CS ratio is also sensitive to the 
gas-phase 0/C ratio. Thus, when observing the SO/CS 
ratio in two individual clouds, it is not easy to disentangle 
the effects of time evolution from those of the 0/C ratio. 



In their survey towards dense cores of dark clouds. 



Suzuki ct al. 1992 found a strong correlation between the 



CCS and HC3N column densities, and, using data from 
Cox et al. 198E, a weaker correlation between CCS and 



C3H2. They explained these correlations by the produc- 
tion of CCS in reactions between hydrocarbons and S^, 
which should be abundant in regions where also carbon- 
chain molecules are formed ( Prasad fc Huntress 1982| ). 
The formation of CS also involves S"*" and occurs in the 
very early stages of chemical evolution. Thereafter, CS 
is mainly recycled via the thioformyl ion, HCS"*", and its 
abundance is considered t o be rather constant in time 
([Nejad fc Wagenblast 199S|) . 



The correlation between HC3N and CCS is due to fact 
that C3H dominates the production of both molecules. 
HC3N is related to C3H2 via the neutral-neutral rcac- 
tion C3H2 + N HC3N + H (jHcrbst fc Leung 1989; 



Nejad fc Wagenblast 1999). On the other hand, complex 
molecules and 



carbon chain 
lation with 



ammonia, 



NH3 



CCS 
(e-g- 



show 



no corre- 



Little et al. 1979 



Suzuki et al. 1992 ). This can be explained by the fact 



The precursor ion of CH3OH, CH3OHJ, is formed 
from the radiative association reaction between CH;^ and 



H2O (e.g. Herbst fc Leimg 198!;). The close relation to 



CH3 makes CH3OH an early-time molecule, as long as 
gas phase reactions are considered. The large abundances 
of CH3OH in high-mass star- forming cores has been 
explained by evaporation from grain surfaces due to 



shock heating (e.g. Menten et al. 1988) 



3. Observations 

The observations were made in February 1998 at the 
Swedish-ESO Submillimetre Telescope (SEST) on La Silla 
in Chile. The 100 and 150 GHz SIS receivers were used 
simultaneously. These were connected to a 2000 chan- 
nel acousto-optical spectrometer (AOS), which was split 
in two bands of 43 MHz each. The resulting channel- 
separations are 0.13 kms~^ (at lOOGHz) and 0.09kms~^ 
(at 150 GHz), respectively. Calibration was achieved by 
the chopper-wheel method. The pointing and focus of the 
telescope were checked at 3-4 hours intervals against the 
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SiO(i' = 1, J = 2 ^ 1) maser line towards AH Sco, U Men 
and R Dor. The pointing accuracy was typically 3" (rms) 
in each coordinate. 

Most of the lines were observed in the frequency 
switching mode with a frequency throw of 6 MHz. Po- 
sition switching was used for CH3CCH and N2H+ since 
the spectrometer band covers several of their lines. 

The observations were started by mapping both 
sources in the HC3N( J =10^9) and S0( Ja, = 43 ^ 82) 
transitions simultaneously in the frequency switching 
mode. HC3N was not detected in CrA C, and therefore the 
map was repeated in the transition pair HNC( J = 1 — > 0) 
& SO(JAf = 43 82). The rest of the observations were 
single position pointings towards the SO (in CrA C) or 
HC3N (in Cha-MMSl) line emission peaks in these maps. 
Table || lists the (0,0) positions used in the mapping 
observations and the peak positions, towards which the 
single pointing observations were made in each core. 

In Table |^ we give a summary of the observed 
lines and some telescope characteristics at the corre- 
sponding frequencies. The columns of this Table are: 
(1) the rest-frequency; (2) the molecule; (8) the tran- 
sition; (4) the line strength of the transition (Sui); (5) 
the energy of the upper level of the transition (Eu/k); 
(6) and (7) the half-power beam-width (6mb) and the 
main beam efficiency (?ymb) of the antenna. Further de- 
tails of the SEST are available e.g. on the W WW at 
http: / /www. Is. eso.org/lasilla/Telescopes/SEST. 



4. Results 

4.1. Maps 

The velocity-integrated HC3N(J = 10 ^ 9) and 
SO{Jn = 43 ^ 82) hne-intensity maps of Cha-MMSl 
are presented in Fig. |^. The locations of the embed- 
ded infrared sources Ced 110 IRS2, IRS4 and IRS6 



( Prusti et al. 1991 ), and the peak p osition of the 1.3 mm 
radio continuum source Cha-MMSl ( Reipurth et al. 1996 ) 
are also indicated. It is evident from the map that the 
HC3N distribution is compact and peaks close to the 
dust continuum maximum, some 30" south of it. A 



comparison with Fig. 1 of Reipurth et al. 1996 shows that 
the dust continuum emission and the HC3N line emission 
have similar sizes and shapes, and therefore conceivably 
originate from the same source, despite the slight shift 
in the peak positions. The infrared point sources lie near 
the northeastern edge of the molecular cloud core. SO 
line emission is present in the compact core, but shows no 
clear concentration. The maxima lie on its northern and 
southern sides, and are possibly related to a molecular 
outflow (see Sect. 5). 

The integrated HNC(J = 1^0) and SO(JAr = 43 ^ 
82) intensity maps of CrA C are shown in Fig. In these 



maps, the situation is the opposite of that in the maps of 
the Cha-MMSl core. HNC, which was chosen instead of 
the undetected HC3N, is weak and dispersed, whereas SO 
has a strong peak towards the offset (80", 30") from the 
centre of the map, and a weaker second maximum towards 
the offset (60", -30"). 

4-2. Spectra and line parameters 

One location in each core was observed in several 
molecular transitions. In Cha-MMSl this position lies 
towards the HC3N peak at the offset (-60", -60"), 17" 
southeast of the dust continuum peak. In CrA C the 
selected position lies close to the SO peak at the offset 
(30", 60") from the map centre, which has the highest SO 
line intensities. The spectra obtained towards these two 
locations are shown in Figs. |^ and ^. 

In Cha-MMSl all observed transitions except 
HC3N(J = 15 — * 14) were detected. The energy of the 
upper state of this transition is rather high (52 K) and 
so the non-detection of this line in a cold cloud is unsur- 
prising. Since the SO lines were rather weak (TJ < 0.7 
K) and showed no self-absorption features, '^'^SO was 
not observed. In CH3CCH, only the = and if = 1 
components of the J = 5 ^ 4 and J = 8 ^ 7 transitions 
were clearly detected. The upper limits (8(7) for the K—2 
and K—3 components are T^=0.063K (J = 5 4) and 
TX=0.045K (J = 8 ^ 7), respectively 

In CrA C, both of the observed SO transitions ( Jat = 
82 2i and 43 82) were bright (TJ > 1.5 K). Only 
a tentative detection of HC3N( J = 10 ^ 9) was made 
at a level of T^=0.05 K, and so thcr transitions were not 
attempted. None of the CH3CCH lines were detected, and 
the rest of the lines, both from carbonaceous species as 
well as from N2H"'", were weaker than in Cha-MMSl. 

The common isotopomer HNC shows a broad, asym- 
metric profile. The fact that the peak is shifted with re- 
spect to the other lines, in particular that of HN^^C, can 
be explained by absorption in the low density halo around 
the core. The larger width of this transition is partly due 
to the hyperfine structure discussed in Sect. 2.1. However, 
the broad wings visible on both sides of the line centre are 
not seen in any other line (including C^^O). The presence 
of high velocity gas towards this location should be tested 
in other optically thick transitions. 

In both clouds the C3H2( Jk^,^-^ = 812 221) spectra 
also contain the Jk = 8_i — ^ 2_i E and Jfe = 3o ^ 2o A+ 
lines of CH3OH. The former component appears in Figs. || 
and H as a negative spike due to the frequency switching 
procedure. 

The observed line parameters (from fits of Gaussian 
line-profiles) in the selected positions are given in Ta- 
ble ^. We list the peak antenna temperatures (TJ), the 
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radial velocities (ulsr)i the full- width at half-maximum 
(FWHM) hue- widths (Aw), and the velocity-integrated 
line-intensities (/ dv) of the detected lines for both 
cores. For the undetected lines, the upper limits (3ct) for 
Tp^ are given. 

4-. 3. Column densities 

The method used in the calculation of the excitation 
temperatures, optical depths, and column densities is 
described in Appendix A. The results of this analysis are 
summarized In Table |[ The columns of this Table are: 
(1) the molecule; (2) the transition; (3), (4) and (5) the 
excitation temperature (Tex) and the line-centre optical 
depth (ti/^, ) of the transition, and the total column 
density (iVtot) of the molecule in Cha-MMSl; (6), (7) 
and (8) Tex, Tj/„, and iVtot for the same molecule and 
transition in CrA C. 

We have estimated the column densities of C^^O to- 
ward s the centres of C rA C and Cha-MMSl using data 
from Harju et al. 1993 and Haikala et al. 2000 (in prepa- 
ration). Assuming optically thin line emission and a uni- 
form excitation of the energy levels with Tex = 10 K, 
we obtained N{C^^O) « 3 10^^ cm'^ in CrA C, and 
N{&^0) « 2.5 10^5 cm-2 in Cha-MMSl. The two sources 
thus appear to have very similar C^^O column densities. 
The column densities of the observed molecules relative 
to C^^O are presented in Table The conversion of the 
ratios listed in this Table to fractional abundances rela- 
tive to H2 is not trivial. The fractional CO abundance is 
expected and has also been observed to change from cloud 



to cloud. As discussed by Harjunpaa & Mattila 1996 the 



fractional C^*0 abundances can be different in star form- 
ing and non-star forming regions. The relations they de- 
rived towards the R CrA and the Cha I clouds (assuming 
N{B.2)/E{J - K) = 5.410^1 cm-2 mag-^) give 2.410^2 
cm-2 and 1.2 10^^ cm'^ for 7V(H2) in CrA C and Cha- 
MMSl, respectively. However, neither region studied by 
Harjunpaa fc Mattila 199^ in these complexes is actually 



close to our sources. The assumption that the relation 
they derived towards Coalsack is valid in the quiescent 
CrA C would result in an H2 column density of 3.8 10^^ 
cm~2. In comparison, the relation between 7V(H2) and 



iV(C^^O) derived by Frcrking et al. 1979| in Taurus yields 



iV(H2) = 1.910^2 cm-2 and iV(H2) = 1.610^2 cm'^ to- 
wards CrA C and Cha-MMSl, respectively. All these H2 
column densities are similar within a factor of two for each 
source, and so the values presented in Table ^ may be con- 
verted to approximate fractional abundances by multiply- 
ing them by 2 X 10^^, bearing in mind that for CrA C, 
the value 1 x 10^^ is possibly more appropriate. 

Table |^ shows that the column densities of carbon 
chain molecules (CCS, HC3N, CH3CCH and C3H2) rel- 
ative to C^^O are much larger towards Cha-MMSl than 
towards CrA C, for which upper limits only could be de- 



rived for aU species other than CCS. The HCgN/Ci^O 
column density ratio towards Cha-MMSl is almost 0.2, 
which implies that the fractional HC3N abundance is of 
order 10~^, being several times larger than towards the 



cyanopolyyne maximum in TMC-1 (Pratap et al. 1997). 
The SO/C^^O abundance ratio is larger towards CrA C. 
Inspection of the intensity ratios of the '^''SO and the nor- 
mal isotope SO lines reveals that the latter are optically 
thick or even self-absorbed towards CrA C, and the de- 
rived SO column densities there are almost certainly lower 
limits (the terrestial '^^S/^^g isotope ratio is 22.5 and in 
the local interstellar medium values larger than 10 are 
normally assumed). On the other hand, the N2H"'"/C^^0 
column density ratio is a few times larger towards Cha- 
MMSl than CrA C. Two molecules which have similar col- 
umn densities towards both cores are C'^^S and CH3OH. 
For the rest of the molecules, the column densities rela- 
tive to C^^O are so different in the two cores that despite 
the uncertainty concerning the iV(H2)/Af(C^*0) conver- 
sion factor, this issue may be discussed below in terms of 
differences in fractional abundances. 

Chemical differences between the two cores are man- 
ifest in Table |^, where we list the column density ratios 
of several molecules. The difference is most marked in the 
HC3N/SO column density ratio, which is at least three or- 
ders of magnitude higher towards Cha-MMSl than CrA C. 
The ratios HC3N/CCS and HC3N/N2H+ are more than 
a hundred times larger, and the C^'*S/"^'*SO ratio is about 
30 times larger towards Cha-MMSl than towards Cra C. 
However, the CCS/N2H"'" column density ratio is similar 
towards both cores. 

5. Discussion 

In the previous Section it became evident that carbon- 
chain molecules are considerably more abundant in Cha- 
MMSl than in CrA C, whereas the opposite is true for 
SO. In particular, the HC3N/SO column density ratio is 
at least 2000 times higher in Cha-MMSl than in CrA C. 
The fact that HC3N, CH3CCH, C-C3H2 and CCS were de- 
tected without difficulty in one of the cores, but were all 
very weak or undetected in the other, is in agreement with 
the common belief that these species thrive in the same re- 
gions. For CCS, the difference in the abundances towards 
the two cores is the smallest for carbon-chain molecules, 
being 'only' about one order of magnitude. 

As discussed in Sect. 2.2, time-dependent models 
of gas-phase chemistry suggest that HC3N and other 
carbon-chain molecules reach their peak abundances at 
early stages of chemical evolution, whereas SO becomes 
abundant first at later times. Adopting this view, the 
HC3N/SO abundance ratios would indicate that Cha- 
MMSl is chemically younger than CrA C. In addition, the 
CS/SO abunda nce ratio has been sugg e sted to depend on 
the cloud age (Bcrgin fc Langer 1997 ; Ruffle et al. 1999 



iNilsson et al 2000|) . Indeed, the CS/SO ratio is clearly 
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larger towards Cha-MMSl, which is mainly due to vari- 
ation of the SO abundance, since the CS abundance is 
rather similar towards both cores. 

On the other hand, it is difficult to understand in terms 
of pure gas-phase chemistry why CrA C has a higher 
SO abundance but a lower N2H+ abundance than Cha- 
MMSl: The production of N2H+ also increases first at 
later stages of chemical evolution. In particular, the fact 
that CCS/N2H+ column density ratio is similar towards 
both cores seems to contradict the expectation presented 
in Sect. 2.2 that this ratio could be usable as a chemical 
clock. 

The inclusion of freezing-out of atoms, ions and 
molecules onto grain surfaces causes a complication for 
the time-dependence of molecular abundances. According 
to Ruffle et al. 1997, a second cyanopolyyne peak may oc- 



time cyanopolyyne peak or a later stage where the chem- 
istry has been influenced by neighbouring young stars. 

The latter conclusion is based on the large N2H+ abun- 
dance in Cha-MMSl. Moreover, the interpretation that 
a collapsing cloud core containing a protostar, as Cha- 
MMSl probably does, has reached the second, 'late-time' 
cyanopolyyne peak seems more natural than that it would 
represent 'early-time' chemistry, since such a core should 



exhibit a high degree of depletion (Rawlings et al. 1992 
Bergin fc Langer 1997| ). Possible outflows and radiation 



from the embedded protostar, IRS 10, or the other nearby 
protostars may have influenced the chemical composition 
of Cha-MMSl. [Suzuki et al. 1992| found that NH3 tends 



to be abundant in star-forming regions, and it seems pos- 
sible that the production of NH3 , and therefore also that 
of N2H+ (see Sect. 2), somehow depend on desorption 
processes. On the other hand, little time has passed since 
the collapse began (which should be coincident with the 
second cyanopolyyne peak) because the duration of the 
In the models of Ruffle et al. 1999, the SO peak comes Class phase is very short (~ 10* years), and it is ques- 
tionable whether the chemistry has really changed signif- 
icantly during this time. 

According to the models of Ruffle et al. 1999| and 



cur when the depletion of gas-phase species becomes sig- 
nificant. The SO abundance also increases considerably 
at this time ( [Bergin fc Langer 1997[ [Ruffle et al. 19991 ). 



slightly after the second HC3N peak (see their Figs. 1 and 
2), probably because the C"^ abundance is then lower. Af- 
ter the SO peak, the dominant species in the gas phase are 
CO, CS and SO, until they too are depleted. It is uncertain 



how comparable the models of Bergin & Langer 1997 and 
Ruffle et al. 1999 are, however, since the former do not 
predict any late cyanopolyyne maximum. 

N2II+ is included in the models of Rawlings et al. 1992 
and Bergin fc Langer 1997] , where depletion onto grains 
are taken into account. In both models N2H"'" shows a 
late peaking, and two reasons for this peaking are given. 
First, N2, which contributes to the N2II+ production, is 
only loosely bound in grain surfaces and returns quickly 
to the gas phase. Second, the drop in the gas-phase abun- 
dances of CO, II2O and other neutral molecules with the 
advance of freezing-out decreases the destruction rate of 
N2H"'", until the electron recombination starts to dominate 
its removal. 

Finally, it should be noted that nearby or embedded 
young stars may revive the ion-molecule chemistry in a 



dense core. As discussed by Hartquist ct al. 1996, an in- 



Bergin fc Langer 1997| , where the chemical development is 
coupled to kinematical models, the compositions observed 
in Cha-MMSl and CrA C correspond to times a few mil- 
lion years after the initial state, which in these studies 
has been defined as a state where hydrogen is molecu- 
lar but all other species are found as atoms or atomic 
ions. If the physical characteristics of the two cores were 
similar, Cha-MMSl would not need to be much younger 
than CrA C. The latter core, however, appears to have a 
lower average density, and therefore the chemical evolution 
has probably been slower there. The result is surprising, 
since star- forming cores of the type Cha-MMSl are nor- 
mally expected to be older than starless cores like CrA C. 
Furthermore, previous observations have shown that the 
carbon-chain molecules are less abundant in star-forming 
regions than in quiescent clouds (Suzuki et al. 1992). 

CrA C has been estimated to be gravitationally 
bound and therefore a potential site of star-formation 



tensified desorption of molecules from grain surfaces, the 
subsequent fragmentation of molecules into atoms in the 
gas-phase, and finally the ionization of atoms in a star- 
forming core may lead to increased abundances of C and 
0+. These could then re-start the synthesis of complex 
molecules, and accordingly the chemistry would have an 
'early-time' character. 

We have seen that chemistry models, regardless of 
whether they take accretion onto grains into account or 
not, predict that SO and N2II+ peak at late stages of 
chemical evolution, and that SO peaks later than IIC3N 
during the pre-stellar phase. Therefore, concerning the 
cores studied here, it is difficult to avoid the following 
conclusions: 1) CrA C has reached an advanced chemical 
stage. 2) Cha-MMSl either represents the so called late- 



(Harju et al. 1993). The present results indicate, however. 



that CrA C has existed long enough to reach an advanced 
chemical stage, which according to chemistry models im- 
plies a minimum age of ~ 10^ years. This longevity raises 
doubts about the star-formation capacity of the core. 

In quiescent cores, star formation occurs due to the 
cloud's self-gravity, which has to work against the ther- 
mal pressure and the pressure and tension of the mag- 
netic fields. Line emission from molecules, especially CO 
(as long as it is abundant in the gas phase), and colli- 
sions between gas and dust particles, can keep a cloud 
nearly isothermal even when contracting, and the im- 
portance of the magnetic support in weakly ionized 
cores has been recognized as the dominant factor pre- 
venting the collapse within the free-fall time, which in 
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the case of a dense core is of the order of 10^ years. 
The magnetic support eventually gives way in the cen- 
tral parts of a core due to ambipolar diffusion, result- 



ing in gravitational instability and collapse (Spitzer 1978 
Mouschovias 197£; Lizano & Shu 1987). The density en- 



hancement required for the gravitational collapse de- 
pends on the initial mass-to-magnetic-flux ratio. The 
rate at which the central density increases is related to 
the timescale of ambipolar diffusion (see e.g. Sect. 2 in 
Ciolek fc Mouschovias 1993 ). The timescale of ambipo- 



lar diffusion in a dense core is proportional to the frac- 
tion al ionization, and is typically of the order o f 10^ years 



(e.g. ppitzer 1978[ [Hartquist fc WiUiams 1989]) . This time 
is shorter than the time needed for reaching a chemi- 
cal equilibrium in pure gas-phase models, but comparable 
with the depletion timescale in models in which accre- 
tion onto grain surfaces are taken into account. Accord- 
ing to Hartquist fc Williams 1989, depletion can influence 



the dynamical state of a core, and in the extreme case, 
even prevent it from collapsing by causing an increase 
in the fractional ionization and removing coolants. Fur- 
thermore, as pointed out by Ciolek fc Mouschovias (1993, 
1994,1995), the properties of dust grains can affect the am- 
bipolar diffusion timescale and thus regulate the evolution 
of a core. 

In simulations by fVanhala fc Cameron 1997 , the state 
of a gas-cloud first becomes non-adiabatic (i.e. the cloud 
is allowed to collapse) when its temperature has risen to 
about 27 K, at which point the rotational levels of H2 be- 
come excited. Even though these authors use perhaps un- 
realistically low CO cooling efficiences (they neglect the 
cooling below 10 K) this model may be qualitatively valid 
for very cold cores with highly reduced abundances of 
heavier molecules. If the temperature is too low, then the 
excitation is insufficient to provide strong enough line- 
emission, and a thermal balance is established. Hence, it 
may be that such cores have a poorer star-formation abil- 
ity than previously supposed. Furthermore, an elevated 
temperature possibly is a useful diagnostic of the onset of 
a collapse, and could be used to identify dense pre-stellar 
cores that are close to forming stars. 

Inferring from its chemical composition, the develop- 
ment in CrA C seem to have been very slow indeed, or 
then the gravitational forces in the core are simply bal- 
anced by thermal and magnetic pressures. In the ambipo- 
lar diffusion model the slow dynamical evolution could be 
explained by 1) a low initial mass-to-magnetic- flux ratio 
in the surrounding cloud; or 2) a high degree of ioniza- 
tion or a high fraction of charged dust grains, which both 
would lengthen the ambipolar diffusion timescale. The fact 
that the N2H+ abundance is lower in CrA C than in Cha- 
MMSl suggests that the fractional ionization is indeed 
higher in the former. An alternative reason for the slow 
progress is that the cooling rate is very low in this core. In- 
dependent determinations of the kinetic temperature, the 
degree of ionization, and the degree of depletion towards 



this core would certainly clarify this issue. Based on the 
results of coupled chemical and kinematical models (e.g. 
Rawlings et al. 1992; Bergin & Langer 1997), the high de- 



gree of depletion inferred from the molecular abundances 
implies a very advanced stage of evolution, and if the core 
is going to collapse, then this has already begun, or will 
do so soon deep in its interior parts. A search for infalling 
motions towards CrA C could be therefore profitable. 
According to Lchtincn et al. 2000| , the region of Cha- 



MMSl may be an example of sequential star- formation, 
in which the core collapse is assisted by external forces 
such as stellar-winds originating from nearby young stars. 
This would conform with a rapid development, and may 
reflect favourable initial conditions for contraction and 
gravitational collapse in the parental cloud. In order to 
judge whether the chemistry in Cha-MMSl is controlled 
by desorption or depletion, it would be useful to observe 
molecules which are predicted to increase in abundance 
when desorption mechanisms are effective. Our tentative 
column density estimates of CH3OH (an example of such 
a molecule) did not, however, show any difference in rela- 
tive abundances between CrA C and Cha-MMSl. Detailed 
chemistry models of the relative timing of the IIC3N, SO 
and N2H+ peaks would therefore be valuable. 

The tiny SO maxima lying symmetrically on both sides 
of the IIC3N core in Cha-MMSl (see Fig. Q) may be re- 
lated to an outflow from the embedded protostar. These 
emission peaks and the Herbig-Haro object HH49/50 lie 
on the same line intersecting the centre of Cha-MMSl, 
and it seems that Cha-MMSl can i ndeed be the driving 
source of HH49/50, as suggested by Rcipurth et al. 1996[ . 
The stronger SO line emission towards the outflow axis 
may be understood in terms of an increased efficiency of 
neutral-neutral reactions involved in the formation of SO 
in shock-heated gas. 

6. Conclusions 

Comparison of observed molecular column densities and 
abundances in two dense cores, Cha-MMSl and CrA C, 
reveal that they have very different chemical compositions, 
which according to time-dependent chemistry models can 
be interpreted as representing different stages of chemical 
evolution. Cha-MMSl, which is dynamically more evolved 



and contains a Class protostar (Rcipurth ct al. 1996 
Lehtinen et al. 2000), has large abundances of carbon- 



chain molecules. This is normally interpreted as a signa- 
ture of an early stage of chemical evolution. In contrast, 
CrA C, which we suspected to be virginal and therefore 
a likely fund of 'early type' molecules, has, to all ap- 
pearances, reached chemical maturity. This result seems 
to contradict the common assumption that star-forming 
cores are (chemically) older than starless cores. The chem- 
ical difference between the two cores is particularly pro- 
nounced in the HC3N/SO abundance ratio, which is about 
2000 times higher in Cha-MMSl than in CrA C. 
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The fact that Cha-MMSl has a large N2H+ abun- 
dance arouses doubts about its chemical youth, however. 
Namely, this molecular ion is formed from N2, which 
should become abundant only at later times. In fact, Cha- 
MMSl probably represents the 'late-time' cyanopolyyne 



peak predicted by Ruffle et al. 1997, which occurs when 



the freezing-out of molecules onto grain surfaces begins to 
be significant. Alternatively, the characteristics of youth- 
ful chemistry in Cha-MMSl could have developed due to 
outflows and ionizing radiation from the embedded pro- 
tostar or the neighbouring newly born stars. It should be 
noted, however, that these two alternatives do not have 
to lie very much apart in time. A high degree of deple- 
tion should be coincident with the collapse of the core nu- 



cleus (e.g. Rawlings et al. 1992 ; Bcrgin fc Langer 1997 ) 



and Cha-MMSl still contains a protostar in its main 
accretion phase. The fact that Cha-MMSl is associated 
with a small cluster of young stars indicates a high star - 
formation efficiency in this region dLehtinen et al. 2000D , 



which would be consistent with a rapid collapse. 

We have detected distinct SO emission maxima on 
both sides of Cha-MMSl, which may indicate the pres- 
ence of a north-south oriented bipolar outfiow. The sug- 



gestion by Rcipurth et al. 1996 that Cha-MMSl contains 



the central source of the Herbig-Haro object IIII49/50 is 
strongly supported by the fact the SO peaks, Cha-MMSl 
and HH49/50 lie along the same line. 

The situation observed in CrA C, i.e. low abun- 
dances of carbon chain molecules and the large SO and 
N2H+ abundances, correspond to final stages of pure gas- 
phase chemistry evolution or to the 'SO peak' which 
comes after the l ate cyanopolyyne peak in the model of 



Ruffle et al. 1999t in which gas-grain interactions has been 
taken into account. The lack of star- formation activity in 
the chemically 'old' core CrA C shows that the dynamical 
evolution has been slow there or that the core may have 
even reached an equilibrium. This may be due to a low 
mass-to-magnetic-fiux ratio in the surrounding cloud, or 
a high degree of ionization or a high fraction of charged 
dust grains making the ambipolar diffusion time very long 



e.g. Ciolek fc Mouschovias 1994). On the other hand, in- 



sufficient cooling efficiency due to freezing-out of heavy 
molecules may have increased the thermal support. 
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Table 1. Spectroscopic properties of the observed molecules. 



Molecule 


Type 


/■^electric 


A 


B 


C 






[Debye] 


[MHz] 


[MHz] 


[MHz] 


SO 


linear (^E) 


1.55 


— 


21 523.020 


— 




linear (^E) 


1.55 


- 


21 102.720 


- 




linear (^E) 


1.957 




24 103.541 




C-C3H2 


asymmetric top (o/p) 


3.43 


35 092.6 


32 212.8 


16 749.1 


CH30H 


asymmetric top (A) 


0.896 


127628.193 


24 687.552 


23 754.594 


CH30H 


asymmetric top (E) 


0.896 


127628.7051 


24 693.0658 


23 757.1371 


CCS 


linear (^E) 


2.81 




6 477.750 




HC3N 


linear (^E,lifs) 


3.724 




4 549.058 




CH3CCH 


symmetric top (A/E) 


0.75 


158 590.0 


8 545.86 


8 545.86 


HNC 


linear (^E,hfs) 
linear (^E,hfs) 
linear (^E,hfs) 


3.05 




45 331.990 




HN"C 


2.699 




43 545.61 




N2H+ 


3.40 




46 586.867 





Table 2. Map centre coordinates and the positions selected for pointed observations. 



Core 


Map centre 


Peak position 




a 5 


Aa 


A<5 




(1950.0) 






Cha-MMSl 


ll''05™28:0 -77°06'32" 


-60" 


-60" 


CrA C 


19''00""33^2 -37°20'22" 


30" 


60" 



Table 3. Observed transitions and some antenna parameters at their frequencies. 



Frcxjiunicy 


Molecule 


Transition 




En/k 




'/mb 


[MHz] 








;k] 


n 




85 338.893 


ortho-c-C3H2 


Jk^,Kc ~ 2i,2 lo,i 


1.500 


4.095'' 


59 


0.75 


85 457.271 


CHaCCH 


Jk = 5o ^ 4o ^ 


5.000 


12.304 


59 


0.75 


85 455.622 


CH3CCH 


Jk = 5i ^ 4i S 


4.800 


11.484*" 


59 


0.75 


87 090.859 


HN^^C 


J{F) = 1(2) ^ 0(1) 


l.OOO'^ 


4.180 


58 


0.75 


90 663.543 


HNC 


J{F) = 1(2) ^ 0(1) 


1.000"= 


4.351 


56 


0.74 


90 978.993 


HC3N 


.7 = 10 ^ 9 


10.000 


24.015 


55 


0.74 


93 173.809 


N2H+ 


J(Fi,F) = 1(2,3) ^0(1,2) 


1.000"= 


4.472 


54 


0.74 


93 870.023 


CCS 


Jn = 87 — > 76 


7.971 


19.892 


54 


0.74 


96 412.940 


c^*s 


J = 2 ^ 1 


2.000 


6.941 


52 


0.73 


97715.388 


3*S0 


Jn = 82 ^ 2i 


2.935 


9.100 


52 


0.73 


99 299.879 


so 


Jn = 32 ^ 2i 


2.933 


9.226 


51 


0.73 


100 076.389 


HC3N 


J= 11 ^ 10 


11.000 


28.818 


50 


0.73 


129 138.898 


so 


Jn = 83 — > 22 


2.670 


25.511 


39 


0.68 


135 775.633 


3^80 


Jn = 4^3 32 


3.935 


15.609 


37 


0.67 


136 464.400 


HC3N 


J = 15 14 


15.000 


52.395 


37 


0.67 


136 728.010 


CH3CCH 


Jk ~ 80 — + 7o A 


8.000 


29.529 


37 


0.67 


136 725.397 


CH3CCH 


Jk = 81 — > 7i i5 


7.975 


28.708'' 


37 


0.67 


138 178.648 


so 


Jn = 43 — > 32 


3.938 


15.857 


36 


0.67 


144 617.109 


c^^s 


J = 3^ 2 


3.000 


13.881 


35 


0.66 


145 089.595 


ortho-c-C3H2 


Jk^,Kc = 3l,2 — > 22,1 


1.250 


13.699^ 


35 


0.66 


145 097.470 


CH3OH 


Jk = 3-1 2-1 E 


2.667 


11.607'" 


35 


0.66 


145 103.230 


CH3OH 


Jk = ^0 2o A'^ 


3.000 


13.929 


35 


0.66 



^ Relative to the ortho ground-state. 

^ Relative to the _E-symmotry species ground-state. 

" This value refers to the J = 1 ^ transition when all hfs components have been accounted for. 
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Table 4. Observational line parameters towards the centres of Clia-MMSl and CrA C. 
The transitions are listed in the same order as in the figures showing the spectra. 



Molecule 


Transition 




Cha-MMSl 






CrA C 












"LSR 


Av 


/ ^idi, 


Tl 


''LSR 


Av 










[K] 


[kms-1] 


[kms-i] 


[Kkms-l] 


[K] 


[km s"-*-] 


[km s"-*-] 


[Kkms^l] 




43 


^32 




_ 


_ 


_ 


0.29±0.02 


5.48±0.01 


0.43±0.03 


0.13±0.02 


SO 


43 


-»32 


0.37±0.04 


4.35±0.02 


0.55±0.04 


0.21±0.04 


1.44±0.04 


5.43±0.01 


0.56±0.01 


0.87±0.04 


so 


33 


^ 22 




_ 


_ 




0.53±0.06 


5.52±0.02 


0.39±0.04 


0.20±0.06 




32 


^ 2i 


_ 


_ 


_ 


_ 


0.50±0.03 


5.51±0.01 


0.57±0.03 


0.33±0.04 


so 


32 


^ 2i 


0.70±0.02 


4.37±0.01 


0.74±0.02 


0.56±0.03 


1.79±0.05 


5.43±0.01 


0.70±0.02 


1.40±0.06 


C34S 


3 - 


-* 2 


O.lOitO.Ol 


4.64±0.02 


0.53±0.04 


0.05±0.01 


0.09±0.02 


5.83±0.02 


0.48±0.05 


0.05±0.01 


CS-'S 


2 -» 1 


0.22±0.01 


4.66±0.01 


0.63±0.03 


0.14±0.01 


0.23±0.01 


5.89±0.01 


0.62±0.03 


0.16±0.02 


CCS 


87 


^76 


0.31±0.02 


4.15±0.01 


0.65±0.03 


0.21±0.02 


0.07±0.01 


5.40±0.03 


0.38±0.06 


0.03±0.02 


0-C-C3H2 


3i 


2 — » 22,1 


0.54±0.06 


4.31±0.02 


0.73±0.04 


0.43±0.05 


<0.26'^ 


_ 


_ 




0-C-C3H2 


2i 


2 l(),i 


1.24±0.02 


4.32±0.01 


0.82±0.01 


1.09±0.02 


0.15±0.01 


5.63±0.01 


0.38±0.03 


0.07±0.01 


CH30H 


3- 


1 ^2_iE 


0.42±0.03 


4.46±0.01 


0.62±0.02 


0.28±0.03 


0.61±0.03 


5.63±0.01 


0.53±0.01 


0.34±0.03 


CH30H 


3o 


^2o A+ 


0.48±0.03 


4.46±0.01 


0.62±0.02 


0.32±0.03 


0.79±0.03 


5.63±0.01 


0.53±0.01 


0.44±0.03 


HN"c 


1 - 













0.09±0.01 


5.70±0.03 


0.76±0.08 


0.07±0.02 


HNC 


1 - 













0.63±0.0^ 


5.21±0.0^' 


o.esdzO.of 


0.57±0.06 


HC3N 


15 


14 


< 0.13'^ 
















HC3N 


11 


10 


0.57±0.03 


4.39±0.01 


0.74±0.03 


0.47±0.04 










HC3N 


10 


9 


0.85 ±0.02 


4.36±0.01 


0.65±0.01 


0.59±0.03 


0.04±0.01 


5.23±0.04 


0.34±0.09 


0.02±0.01 


CH3CCH 


80 


^7o 


0.20±0.02 


4.53±0.02 


0.58±0.04 


0.12±0.02 


<0.083'' 








CH3CCH 


81 


^ 7i 


0.16±0.02 


4.49±0.02 


0.52±0.04 


0.08±0.02 


<0.083'' 








CH3CCH 


5o 


^4o 


0.43±0.01 


4.41±0.01 


0.63±0.02 


0.30±0.03 


<0.071* 








CH3CCH 


5i 


^4i 


0.39±0.01 


4.41±0.01 


0.62±0.02 


0.26±0.02 


<0.071* 








N2H+ 


1 - 





1.67±0.03 


4.57±0.01 


0.63±0.01 


5.02±0.03 


0.45±0.02 


5.77±0.01 


0.46±0.01 


0.94±0.03 



''30-. 



Applies for the strongest component in the spectrum. 

Table 5. The excitation temperatures and the line-centre optical depths of the observed transitions, and the total 

molecular column densities towards the peak positions of Cha-MMSl and CrA C. 



Molecule 


Transition 




Cha-MMSl 






CrA C 










Tex 




Ntat 








A^tot 








[K] 




[cm-' 


] 


[K] 




[cm"' 


] 


34SO 


43 ^ 32 










3.9±0.1 


0.69±0.01 


1.9±0.1 • 


1013 


SO 


43 ^32 


4.4±0.1 


0.59±0.07 


2.0±0.2 • 


1013 


5.8±0.9 


2.10±0.70 


6.6±2.3 • 


1013 


SO 


33-^22 










5.8±0.9 


0.37±0.05 


7.2±2.2 ■ 


1013 


34SO 


32 ^2i 










3.9±0.1 


1.16±0.01 


2.0±0.1 ■ 


1013 


SO 


32 -^2i 


4.4±0.1 


1.10±0.20 


2.7±0.3 • 


1013 


5.8±0.9 


2.30±0.90 


6.7±2.7 ■ 


1013 


C34S 


3 2 


3.3±0.1 


0.49±0.04 


5.7±0.5 • 


IOI2 


3.4±0.1 


0.38±0.04 


3.9±0.5 • 


IOI2 


C^^S 


2^1 


3.3±0.1 


l.OOitO.lO 


5.7±0.6 • 


IOI2 


3.4±0.1 


0.83±0.12 


5.0±0.7 • 


IOI2 


CCS 


87^76 


4-l-ia 




3.1±0.3 • 


10« 


4±ia 




4.0±2.0 ■ 


IOI2 


0-C-C3H2 


3l,2 22,1 


4.7±0.2 


0.83±0.14 


4.7±0.8 • 


1013 










0-C-C3H2 


2l,2 —> lo,l 


4.7±0.2 


2.80±1.60 


4.7±2.7 • 


1013 


4±ia 




1.2±0.2 ■ 


IOI2 


CH30H 


3^1 2-1 E 


4±ia 




6.8±1.3 • 


1013 


4±ia 




8.0±2.0 • 


1013 


CH30H 


3o ^ 2(, A+ 


4±ia 




1.2±0.2 • 


1014 


4±ia 




1.7±0.3 • 


10" 


HNi^c 


1^0 










4-l-ia 




3.2±0.7 • 


toll 


HNC 


1^0 










4±ia 




2.0±0.1 • 


1012 


HC3N 


11 ->■ 10 


4.1±0.1 


1.10±0.20 


4.5±0.8 • 


1014 










HC3N 


10 ^ 9 


4.1±0.1 


3.00±2.00 


4.5±2.5 • 


10" 


6±1'' 




< 5.7 • lOii b 


CH3CCH 


80 7o 


9.9±0.8 


0.05±0.01 


1.1±0.2 • 


10" 










CH3CCH 


81 ^ 7i 


8.5±0.8 


0.05±0.01 


1.0±0.2 ■ 


10" 










CH3CCH 


5o 4o 


9.9±0.8 


0.09±0.01 


1.1±0.2 • 


10" 










CH3CCH 


5i ^4i 


8.5±0.8 


O.lOitO.Ol 


1.0±0.2 • 


10" 










N2H+ 


1^0 


4±1^ 




1.4±0.1 • 


1013 


4±ia 




2.6±0.1 • 


1012 



* The excitation temperature has been assumed. 



The rms is so large compared to the line intensity that only a upper limit for the columndensity is given here. 
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Table 6. Column densities relative to C^^O*. 



Molecules 


Cha-MMSl 


CrA C 




A^tot/iVci«o (xlO^) 


iVtot/iVciBo (xlO') 




- 


6.3 ±0.5 


SO 


8.0 ± 1.0 


22 ±8 




2.3 ±0.3 


1.3 ±0.2 


CCS 


12 ±2 


1.3 ±0.7 


0-C-C3H2 


19 ±4 


< 0.4 


CH30H 


27 ±6 


27 ±7 


HN"C 




0.11 ±0.02 


HNC 




0.67 ±0.06 


HC3N 


180 ± 35 


< 0.2 


CH3CCH 


44 ±9 




N2H+ 


5.6 ±0.6 


0.87 ±0.07 



* The C^*0/H2 column density ratio is normally assumed to be 
about 2-10~^. In quiescent regions like CrA C this ratio can be , 
however, lower by a factor of two ( Harjunpaa fc Mattila 199E| ). 



Table 7. Column density ratios. 



Molecules 


Cha-MMSl 


CrA C 


HC3N/CCS 


15 ±3 


< 0.14 


HC3N/SO 


23 ±5 


< 0.01 


HC3N/N2H+ 


32 ±6 


< 0.22 


HC3N/CH3CCH 


4± 1 




HC3N/C3H2 


10 ±3 


< 0.48 


C^^S / 3*SO 


6.4 ±0.9=^ 


0.21 ±0.03 


N2H+/S0 


0.7 ±0.1 


0.04 ±0.01 


CCS/N2H+ 


2.2 ±0.3 


1.5 ±0.8 


C=*'*S/N2H+ 


0.41 ±0.05 


1.5 ±0.2 



The ^^S/^'^S isotope ratio is assumed to be 22.5. 
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Cha-MMSl HC3N(10^9) 3.5-5.5 km/s 



- 



-1 - 



2 








-3 



Cha-MMSl S0(4o^3p) 3.5-5.5 km/s 







<] 



-2 








-lAa[']-2 -3 



Fig. 1. The HC3N(J 10 ^ 9) and SO(Jjv = 43 82) integrated intensity (/ Tldv) maps of Cha-MMSl. The 
velocity range and the intensity levels (in Kkms^^) are indicated in the Figure. Typical 3a rms noise levels are 0.04 
Kkrns"^ and O.llKkms"^ in the HC3N and the SO maps, respectively. The locations of the 1.3 mm dust continuum 
peak MMSl (Reipurth et al. 1996), and the nearest embedded infrared sources (Prusti et al. 1991) are denoted by a 
cross and stars, respectively. The beamsizes are indicated in the bottom right of each map. The coordinates of the 
(0, 0) position are given in Table |^. 
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CrA C HNC(1 ^ 0) 4.5-6,5 km/s 



2 



1 - 



<i 



- 




3 



2 



OAa['] -1 



-2 



CrA C 30(43^32) 4.5-6.5 km/s 



2 



1 A 



<i 







1 - 




3 



2 



OAa['] -1 



-2 



Fig. 2. The HNC(J = 1^0) and SO(Jjv = 43 82) integrated intensity (/ T^dv) maps of CrA C. The velocity 
range and the intensity levels (in Kkms~^) are indicated in the Figure. Typical Str rms noise levels are 0.06 Kkms"^ 
and O.llKkms""'^ in the HNC and the SO maps, respectively. The beamsizes are indicated in the bottom right of each 
map, and the coordinates of the (0, 0) position are given in Table |. 
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Cha-MMSl 



3 



2 









CCS ^ 


1 (Jjj=87^76) 


"SO n (Jn=43-32) - 




\ (JK„K„=3jg->2gi) 


SO ^ (Jn=32^2i) 




(Jk^,k„^2i2^ Iqi) 




HCgN 


(J=15^14) 


-C^*S (J=3^2) 


HC3N 


\ (J=11^10)- 


C^*S „ (J=2^1) 

x2 r \ 


HC3N 


1 (J=10^9) 



0246802468 10 
1.5 I 



0. 




-5 5 10 -20 20 40 



Velocity [km/s] 

Fig. 3. The observed spectra towards the centre of Cha-MMSl. The spectra are in the scale and the velocity in 
relative to the local standard of rest. The upper panel shows the lines observed using frequency switching, while the 
lower panel shows the position switched spectra. Note the different velocity ranges. The intensity in the C^^S spectra 
are multiplied by a factor of two. The feature showing a negative intensity in the C3H2 (Jko.aTc = 3i,2 '^2,1) spectrum 
is caused by the Jfe = 3o ^ 2qA'^ line of CH3OH, which lies close in the frequency. The frequency switching method 
causes the line to appear with a negative intensity. The CH3CCH spectra are centered on the K = 2 components, 
which, however, remain undetected. The locations of the K = 0,1,2 and 3 components are marked by arrows. 
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CrA C 



1.5 



0.5 



x2 




2 4 6 



2 4 6 8 10 



(J=1^0) 




CHgCCH 

K= 1 



(J = 8^7) 



CHgCCH (J=5^4) 



i i 



15 



-20 



20 



40 



Velocity [km/s] 



Fig. 4. The observed spectra towards the centre of CrA C. The spectra arc in the scale and the velocity in relative 
to the local standard of rest. As for Cha-MMSl, the frequency switched and position swithed spectra are placed into 
different boxes. The selection of transitions is, however, slightly different. Note that many of the spectra are either 
multiplied or divided by a factor of two. The C3H2( Jk^^k^ = 3i2 —>■ 22i) spectrum contains the negative 'CH3OH 
feature'. In the CH3CCH spectra we have marked the positions where the K = 0,1,2 and 3 components should lie. 
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Appendix A: Formulae used in the excitation 
temperature and column density 
calculations 

The excitation temperatures (Tex) and molecular column 
densities (-/Vtot) were estimated using a local thermody- 
namic equilibrium (LTE) approach in the sense that we 
assume a uniformly excited and Boltzmann distributed 
population of the energy levels for the molecule question. 
However, we refrain from making any assumptions 
concerning the (line-centre) optical depths {t^^i ) . Instead 
we derive an exact equation (A5 below) for a two-level 
system, from which Tex can be determined numerically 



(see Nummclin ct al. 200C for an LTE treatment of the 
excitation of molecules with numerous observed lines 
with no a priori assumptions about the optical depths). 

By assuming that Tj^^j has a gaussian velocity distribu- 
tion, the formal definition of the optical depth leads to 
the following expression for its value at the line-centre: 



Eqs.(A2 and A3) yield: 



In 



[Fivuh Tex) - Fjv^i, Tbg)] 

kTB 



[F(l/ul,Te,) -F(l/ui,Tbg)] 



(A4) 



Combining the Eqs. (Al and A4) and re-arranging the 
terms leads to an expression, which besides Tex contains 
only known constants or observed quantities. By applying 
this to two spectral lines of the same molecule, and taking 
the ratio of the resulting expressions, the final equation 



In 



[f(-^2.^ex)--F(-^2.^bg)] 



[F{.x2,Te,x)-F(.^2,T'bg)]- 



[F(vi,Tcx)--F(>'l,Tbg)J 



F(^l,Tcx)--F'(>'l,rbg) 



F(^2,Tex) 



F(iyi,Tcx) 



exp(- 



A^;iS„i_2 

Al>2Sul,l 



(A5) 



10 



-41 



In 2 167r3 5ui.g(/, X)iVtot 



TT 3h 

X exp(-£;u/kTcx; 



Q{Tex)Av 

(exp(hz^ui/kTc; 



1) 



(Al) 



The factor 10~^^ stems from using Debye and kms~^ as 
the units for the electric dipole moment (fJ-ei.) and the 
FWHM line- width (Au), respectively. Quantities related 
to the observed transition are: the rest-frequency v^i [Hz], 
the line-strength S^i, and the energy of the upper level 
-Eu. Furthermore, A'tot is the total column density [cm"^], 
(5(Tex) is the partition function, g(I, K) is the statistical 
weight due to possible nuclear spin and iC-degeneracy, 
and h and k are the Planck and Boltzmann constants 
(given in cgs- units), respectively. 

A second expression for t^^^ is obtained from the "antenna 
equation" : 



The subscripts "1" and "2" refer to "hne 1" and "line 2" , 
respectively. Terms involving Tex have been collected on 
the left hand side of this expression. We solved Eq.(A5) 
for given observed line parameters by searching through 
a sufficiently large Tex interval. When applying this 
equation to the observed data, we found that in addition 
to the case with a single unique solution, there are also 
cases where no solution exists, as well as cases where 
two solutions are present. The lack of solutions indicates 
that the assumption of a uniform Tex for both transitions 
has been violated or that the beam-filling factors of the 
two transitions are different; in this case we have chosen 
the point where we have the least deviation from a 
solution as the value for Tex. The case with two solutions 
manifestates itself by giving one solution with a low Tex 
and high Ti^^j, the second solution having a higher Tex 
and lower r^^j . In the analysis we have chosen to use the 
low Tex - high Tjy^^j alternative. 



Tr = 



hi^ui 



ri?/ \ T?i M ri „ ( M/'Ao^ Once Tex was estimated, , was calculated from Eq.(A4) 

[i* (I^ul, Jexj - -f' (.I'ub-tbgjJ [1 - exp(-ry^JJ(A^j ^ • jr -c^ 



where 



T(z.ui,T) 



exp (hi^ui/kT) — 1 



(A3) 



Tb is the brightness temperature calculated from the ob- 
served T^ by correcting for relevant beam-efficiencies and 
beam-filling factors. In the analysis of the observational 
data, we used the main-beam brightness temperature 
(Tmb=T^/'7mb), assumcd the beam-filling to be one for 
all lines, and did not introduce any corrections due 
to different beam sizes.. Tbg is the temperature of the 
background radiation (here taken as a black-body at a 
temperature of 2.728 K). 



and subsequently TVtot was derived from Eq.(Al). 

In cases where we have only one observed transition, we 
have assigned a Tex close to that of a species with similar 
excitation requirements. N^o\ is then obtained, assuming 
optically thin emission, from the formula: 



iVt, 



1.67 • 10^3 exp ( Eu/kTex)Q(Tex)/li, 

viilS^,g{I,K) 



1 



(A6) 



Here /iinc= / Tb dv is the velocity- integrated line- intensity 
[K kms-i]. 

In the calculation of the column densities the 
partition functions need to be evaluated. For lin- 
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ear rotors (^T,) we used the following formula (e.g. 
Townes & Schawlow 1975, chapter 1): 



Q 



hB 3 
where B is the rotation constant. 



(A7) 



For C-C3H2 we used: 



(A13) 



h3 ABC 

The contributions from both ortho- and para-states have 
been accounted for by taking into account their different 
statistical weigths: 3 for ortho and 1 for para. 



For CCS and SO we made a direct summation over the 
energy levels with Ei/k < 575 K: 



Q = ^ ffi exp (-E'i/kTcx) 



(A8) 



where gi is the statistical weigth, 2Ji + 1, with the 
relevant J, used. 



For CH3OH we used: 



(A14) 



h3 ABC 

The (equal) contributions from both the A- and the 
_B-symmetry species have been accounted for. 

For CH3CCH it is common to use the expression: 



We would like to emphasize that the evaluation of the 
partition function through direct summation is essential 
in the case of SO and CCS. For example for SO, the 
"usual" integral approximation of the partition function 
(3 — 1) overestimates Q by a factor of 1.1 to 1.85 
in the temperature range 2-70 K; the largest discrepancy 
occurs at Tex ~ 5K. Thus for SO this simple approxima- 
tion can be used only for temperatures above 80 K or so. 
A better analytical approximation for Q at temperatures 
below 80 K can be obtained by treating each of the three 
N ladders as separate rigid rotors and integrating over J. 
Note that only one of the ladders ( J = — 1) has J=0 
as its lowest level; the two other ladders begin with J=l. 
The resulting partition functions are: 



Qj= 



iV-l. ladder 



/J=Af,laddcr 



^J=Ar+l, ladder 



k Tex 

hB 



(A9) 



^ - ^ ) exp(-A£;i/fcrex) (AlO) 



^-1) eM-^E^/kT,^) (All) 



where AEi is the energy difference between the lowest 
levels in the J — N and J = A^ — 1 ladders, and AE2 
the corresponding difference between the J = A^ -I- 1 and 
J = N — 1 ladders. The lowest states in these ladders are 
Jn=Oi (the ground-state), li, and Iq, respectively. 

The total partition fimction is then: 

Q — Q,J=N~1, ladder + Qj=JV, ladder + Qj=JV+l, ladder • (A12) 

For SO, this approximation overestimates the true parti- 
tion function by less than 15% for temperatures between 
3 and 80 K, the error being less than 5% for Tex > 25 K. 
The values for AEi and AE2 are 15.18K and 1.44K 
(SO); 9.33 K and 0.53 K (CCS), respectively. 



8 Uk^Tl 



S^h^ABC-'-'''^-- 
where the contributions from both the A- and the 
i?-symmetry species have been accounted for by taking 
into consideration the spin statistics of the various levels 



(see e.g. Townes & Schawlow 1975, chapter 3) 



However, Eq.(Al5) is valid only at high enough tem- 
peratures. Since the lines observed by have low excitation 
temperatures, we derived a more accurate approximation 
for temperatures below 30 K. Each if-ladder was treated 
as a linear rotor and the integration was performed over J. 
Note that J=K is the lowest level of a specific AT-ladder. 
By taking into accoimt only the three first ladders {K—0, 
1, 2) we have: 

k Tex, J 4 

3 



;K=0, ladder 



Qk= 



Qk- 



1, ladder ' 



:2. ladder 



hB 
hB 



(A16) 



- ) exp(-A^3/fcTex,if) (A17) 



- ^ - 12 exp(-A£;3/A;Tex„7) ) x 



X exp(-A£;4//cTex^K) 



(A18) 



where A£;3/k « 7.20 K and AE^/k « 28.80 K are the 
energy differences between the lowest levels of the K—1 
and K—2 ladders, respectively, relative to the A'=0 lad- 
der. Tox.j is the excitation temperature along a specific 
AT-ladder (i.e. K fix, J varies), and Tcx,k is the excitation 
temperature between different A'-ladders (i.e. K varies, J 
fix). 

The resulting expression for the total partition func- 
tion for CH3CCH, which in the indicated temperature 
range underestimates the result from a direct summation 
by less than 10%, is given by: 



ladder 



ladder 



^K^2, ladder 



(A19) 
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